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Abstract. The cosmic far-infrared background is now well measured 
| from 140 /im to 1 mm. Uncertainties remain at 100 fim (and even more 

. at 60 /im). These are dominated by limitations of the zodiacal model. 

| The nature of sources dominating the background near its maximum are 

beginning to be studied through deep surveys carried out with ISOPHOT. 
00 | These surveys show very steep number counts and fluctuations of the 

background due to unresolved sources. 
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1. The Cosmic Far-InfraRed Background (CFIRB) 



1.1. Determination of the CFIRB above 100 /ira: another Approach 

In the parts of the sky with no molecular clouds or dense HII regions, the far-IR 
emission can be written as the sum of dust emission associated with the neutral 
Q-f gas and with the diffuse ionised gas, interplanetary dust emission, the CFIRB, 

and the cosmic microwave background and its dipole. In previous studies (Puget 
et al. 1996; Fixsen et al. 1998; Hauser et al. 1998), dust emission associated 
with the ionised gas which could not be traced independently has been either 
not subtracted properly or neglected. 

Lagache et al. (1999) have for the first time detected dust emission in the 
ionised gas and shown that the emissivity (which is the IR emission normalised 
to unit hydrogen column density) of dust in the ionised gas is nearly the same as 
in the neutral gas. This has consequences for the determination of the CFIRB. 
Following this first detection, Lagache et al. (2000) have combined HI and 
WHAM H Q data (Reynolds et al. 1998) with far-IR COBE data in order to 
derive dust properties in the diffuse ionised gas as well as to make a proper 
determination of the CFIRB. After careful pixel selection they describe the far- 
infrared dust emission as a function of the HI and H + column densities by: 

IR = A x N(HI) 2Qcm -2 +Bx N(H+) 2Qcm -2 + C (1) 

where N(HI) 2 o C m- 2 an d N(H + ) 2 o C m- 2 are the column densities normalised to 10 20 
hydrogen atoms and ions per square centimeter respectively. The coefficients A, 
B and the constant term C are determined simultaneously using regression fits. 
They show that about 25% of the IR emission comes from dust associated with 
the diffuse ionised gas which is in very good agreement with the first deter- 
mination of Lagache et al. (1999). The CFIRB spectrum obtained using this 
galactic far-infrared emission decomposition is shown in Figure 1 together with 
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Figure 1. CFIRB spectra obtained from the decomposition of the far- 
infrared sky (continuous line) and determined for the Lockman Hole 
region (dashed line) by Lagache et al. (1999). Also reported are DIRBE 
values at 100, 140 and 240 fim from Lagache et al. (2000). 



the CFIRB FIRAS determination of Lagache et al. (1999) in the Lockman Hole 
region. We see very good agreement between the two spectra, as well with the 
Fixsen et al. (1998) determination (Figure 2). Fixsen et al. (1998) used an in- 
frared tracer for the interstellar component which is equivalent as long as the 
dust emission spectrum from the diffuse neutral and ionised interstellar gas are 
very similar. This hypothesis is confirmed by the Lagache et al. (2000) analysis. 

At 140 and 240 ^m, the values obtained for the CFIRB are 1.13±0.54 MJy/sr 
and 0.88±0.55 MJy/sr respectively. For each selected pixel, Lagache et al. (2000) 
compute the residual emission, R = IR - AxN(HI) - BxN(H + ). Uncertainties 
in the CFIRB have been derived from the width of the histogram of R (sta- 
tistical uncertainties derived from the regression analysis are negligible). The 
CFIRB values obtained, although much more noisy (due to the small fraction of 
the sky used), are in very good agreement with the determination of Hauser et 
al. (1998). At 140 /im, the CFIRB value of Lagache et al. (1999) is smaller than 
that derived here since the assumed WIM (Warm Ionised gas) dust spectrum 
was overestimated (the WIM dust spectrum was very noisy below 200 /im and 
the estimated dust temperature was higher). At 100 ^m, assuming an accurate 
subtraction of the zodiacal emission, the far-IR emission decomposition gives: 
Icf.trb(100)= 0.78±0.21 MJy/sr. This is the first time that two independent 
gas tracers for the HI and the H + have been used to determine the background at 
100 fim. The CFIRB value of 0.78 MJy/sr can be compared to the non-isotropic 
residual emission found by Hauser et al. (1998). The average over three regions 
of the residual emission, equal to 0.73±0.20 MJy/sr, is in very good agreement 
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Figure 2. Comparison of the CFIRB determinations: Lagache et al. 
2000 (jagged continuous line and light diamonds), Hauser et al. 1998 
(dark diamonds, arbitrarily shifted), Fixsen et al. 1998 (dotted and 
smooth continuous lines) 

with the determination of Lagache et al. (2000). 

So we see, using different approaches, that we are now converging on the 
shape and level of the CFIRB above 100 jim. 

1.2. Where are the Main Uncertainties? 

There are two main uncertainties in the CIRB (Cosmic InfraRed Background) 
determination in the IR and sub-mm domain: (1) Dust emission associated with 
the ionised gas and (2) Zodiacal emission. 

If the spectrum of dust associated with the ionised and neutral gas is the 
same, a CFIRB determination using color ratios should give the same values as 
other methods implying independent HI and H Q gas tracers. However, one has 
to note that methods based on the intercept of the far-IR/HI correlation for the 
determination of the CFIRB are dangerous. For example, for the parts of the 
sky selected by Lagache et al. (2000), this intercept is about 0.91 MJy/sr, which 
is quite different from the value of the CFIRB. 

The zodiacal emission was obtained by Kelsall et al. (1998) relying on its 
time variability. Other models have been built for the zodiacal emission which 
is critical to a determination of the CIRB in the near infrared (Wright 2000; 
Finkbeiner et al. 2000). The accuracy of these models can be estimated by 
comparing the residuals observed at wavelengths were the zodiacal emission 
is maximum (12 and 25 fj,m). For example, on the one hand, the Kelsall et 
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al. (1998) model leads to residuals of about 4.7 x 10~ 7 Wm _2 sr _1 at these 
wavelengths. On the other hand, upper limits on the CIRB have been established 
in this wavelength range using TeV gamma rays from extragalactic sources at 
the level of 10~ 8 Wm~ 2 sr _1 (e.g., Biller et al. 1995). This is not in contradiction 
with the estimate of uncertainties in the Kelsall et al. zodiacal model by Hauser 
et al. (1998), who give only upper limits for the CIRB. Nevertheless this result, 
together with the fact that the residuals follow a zodiacal emission spectrum, 
suggests that this best model underestimates the zodiacal emission as the other 
contributions (instrumental and interstellar) are significantly smaller. We can 
thus make a conservative estimate of the amount of zodiacal emission which 
is not removed by this model at 12 and 25 /im to be 4 x 10 -7 Wm _2 sr _1 . 
The amounts not removed at 60 and 100 /im are thus 4 x 10 -8 Wm _2 sr _1 and 
8.4 x 10 -9 Wm~ 2 sr _1 (using the Kelsall et al. smooth high latitude zodiacal 
cloud color ratios). This reduces the CIRB at 100 fim from 0.78 to 0.50 MJy/sr. 
At 60 //m the extra zodiacal emission to be removed (4 x 10~ 8 Wm _2 sr _1 ) is 
comparable to the residuals, and thus no meaningful value can be obtained on 
the CIRB at this wavelength. 



2. Resolving the CFIRB at 170 /xm 
2.1. The FIRBACK Deep Survey 

FIRBACK is a survey of 4 square degrees in 3 high galactic latitude fields, cho- 
sen to have as low a HI column-density as possible, typically Nhi — 10 20 cm~ 2 , 
and as much multiwavelength coverage as possible. Observations were carried 
out with the ESA Infrared Space Observatory, ISO, (Kessler et al. 1996) us- 
ing the ISOPHOT photometer (Lemke et al. 1996) with the C200 camera and 
Ci6o broadband filter centered at A = 170 ^m. A detailed description of the 
reduction, data processing, and calibration will be discussed in Lagache & Dole 
(2001); the analysis of the complete survey is discussed in Dole et al. (2001). 

Eighty-six sources are detected above the sensitivity limit of the survey 
(200 mJy, 5a). The number of sources detected above 120 mJy (3<r) is 235. The 
first result of this survey is the high number of sources observed when compared 
to no, or moderate, evolution models for infrared galaxies. Extensive simulations 
were done to establish noise properties, incompleteness (80% at 200 mJy) and 
Eddington bias; the cumulative number counts were then established from this 
catalog. These number counts are extremely steep, N(> S) oc S~ 2 ' 2 , indicating 
very strong cosmological evolution of infrared galaxies. This is illustrated in Fig- 
ure 3(taken from Dole et al. 2000). Two models from Guiderdoni et al. (1998) 
are shown. It can be noticed that the semi-empirical model "E" , which was built 
to explain the infrared background, fits the counts well but falls slightly short 
of explaining the steepness of the counts. 

Number counts in a small field (0.25 square degrees) have already been 
published (Puget et al. 1999). This field was observed as a "feasibility demon- 
stration" of the FIRBACK survey at a time when the limitations to weak source 
detection with ISOPHOT were not yet well understood. The counts given in 
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Figure 3. FIRBACK galaxy counts superposed on the Guiderdoni 
et al. (1998) models "A" (dotted line) and "E" (continuous line) (from 
Dole et al. 2000). 



Puget et al. (1999) are compatible with the counts shown in Fig. 3. 

Integrating these counts for all sources detected above 120 mJy (3a) gives a 
brightness of 1.6 nW raT 2 sr _1 or 8% of the CIRB at 170 /j,m. If we extrapolate 
the counts to weaker fluxes with a slope -2, we find that the background is fully 
accounted for by sources brighter than about 10 mJy. As we expect the counts 
to flatten in a progressive way we predict that the background at 170 //m is 
likely to be dominated by sources of a few mJy. 

Considering the rather low angular resolution allowed by a 60 cm telescope 
at this wavelength, the identification of these sources with optically detected 
sources is difficult. Starburst galaxies are expected to be radio sources with a 
well defined ratio of infrared to radio flux. Radio surveys have been carried 
out in the FIRBACK fields with the VLA for the two northern fields (Ciliegi 
et al. 1998) and the Australian Telescope for the southern field. Furthermore, 
15 jj,vn observations have been carried out with the ISOCAM instrument aboard 
ISO (Oliver et al. 2000). These two sets of observations are the best tools for 
initiating the identification process when counterparts are available — which is 
the case for about 50% of our sources. Confirming the identification is best done 
with ground based millimeter or submillimeter observations with SCUBA, CSO 
or the IRAM 30 meter telescope. The last step is optical identification using 
the radio positions if a millimeter or submillimeter detection has confirmed the 
FIRBACK-radio tentative identification. This process is frustratingly slow. On 
the one hand, a substantial fraction of the sources are relatively nearby moderate 
starburst galaxies with bright optical counterparts (redshift less than 0.3, far 
infrared luminosities of a few 10 11 L Q ). On the other hand, a small fraction 
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of the sources have very weak or no optical counterpart and are likely to be 
associated with distant (redshift larger than 1) ultraluminous starburst or dust 
enshrouded AGNs. 

2.2. The Lockman Hole Survey 

A survey of comparable depth has been carried out in the Lockman hole (Kawara 
et al. 1998; Matsuhara et al. 2000). The number counts at 170 /im are in excellent 
agreement with the counts from the FIRBACK survey. It should be noticed that 
the Lockman hole is one of the fields with the lowest cirrus content. The agree- 
ment in the number counts confirms that, for point sources, the contamination 
of the extragalactic source counts by small scale cirrus structure is negligible. 
This is also in agreement with the conclusions of the analysis reported in the 
next section, which discusses the fluctuations of the background. The difference 
in the power spectra of the cirrus component and of the extragalactic component 
leads to a negligible cirrus contribution to the power spectrum at the highest 
spatial frequencies explored by this survey. 

3. The 170 CFIRB Fluctuations 

3.1. Why Search for the CFIRB Fluctuations? 

The CFIRB is taken to consist of sources with number counts as a function of 
flux which can be represented, for the present discussion, by a simple power law: 



Obviously, these number counts need to flatten at low fluxes to insure a finite 
value of the background. Thus, we assume that a=0 for S < S*, where S* is 
the flux where the counts flatten. 

For the simple euclidian case (a=1.5), the CFIRB integral is dominated 
by sources near S* and its fluctuations are dominated by sources which are 
just below the detection limit Sq. It is well known that strong cosmological 
evolution, associated with a strong negative K-correction, could lead to a very 
steep number count distribution (see, for example, Guiderdoni et al. 1998 and 
Franceschini et al. 1998). The present far-IR observations indeed show a very 
steep slope of a=2.2 (Dole et al. 2000). In this case, the CFIRB integral is still 
dominated by sources near S*, but its fluctuations are now also dominated by 
sources close to S*. Thus, it is essential to study the extragalactic background 
fluctuations, which are likely to have a substantial contribution from sources 
with a flux comparable to those dominating the CFIRB intensity (sources that 
are not just below the sensitivity limit). 

To see if we can detect the CFIRB fluctuations, we need wide field far-IR 
observations with high angular resolution and very high signal to noise ratio. 
The FIRBACK project fulfills all these conditions. To search for CFIRB fluctu- 
ations, we have first used the so-called "Marano 1" field, for which we obtained 
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4 independent coadded maps which allow us to properly determine the instru- 
mental noise (Lagache & Puget 2000). In this field, we have a signal to noise 
ratio of about 300 and we detect 24 sources that we remove from the original 
map. We first carry out an analysis on this field, and then extend it to the other 
FIRBACK fields. 



3.2. CFIRB Fluctuation Mean Level 

FIRBACK maps show background fluctuations which consist of two components 
that we want to separate, galactic cirrus fluctuations and, if present, extragalac- 
tic fluctuations. 

Our separation of the extragalactic and galactic fluctuations is based on a 
power spectrum decomposition. This method allows us to discriminate between 
the two components using the statistical properties of their spatial behaviour. 
Fig. 4 shows the power spectrum of the "Marano 1" field. In the plane of the 
detector, the power spectrum measured on the map can be expressed in the 
form: 

Pmap = Pnoise ~i~ {Pdrrus ~\~ Psources) X Wfc (3) 

where Pnoise is the instrumental noise power spectrum measured using the 4 
independent maps of the Marano 1 field (Lagache & Puget 2000), Pdrrus anci 
Psources are the cirrus and unresolved extragalactic source power spectra respec- 
tively, and Wk is the power spectrum of the point spread function, PSF. For our 
analysis, we remove Pnoise from P m ap- 

We know from previous work that the cirrus far-infrared emission power 
spectrum has a steep slope, Pdrrus (Gautier et al. 1992; Kogut et al. 1996; 

Herbstmeier et al. 1998; Wright 1998; Schlegel et al. 1998; Miville-Deschenes et 
al. in preparation). These observations cover the relevant spatial frequency range 
and have been recently extended up to 1 arcmin using HI data of very diffuse 
regions (Miville-Deschenes 1999). The extragalactic component is unknown but 
certainly much flatter (see the discussion in Lagache & Puget 2000). We thus 
conclude that the steep spectrum observed in our data at k <0.15 arcmin -1 
(Fig. 4) can only be due to cirrus emission. The break in the power spectrum 
at k ~ 0.2 arcmin -1 is very unlikely to be due to the cirrus emission itself which 
is known not to exhibit any preferred scale (Falgarone et al. 1998). Thus, the 
normalisation of our cirrus power spectrum Pdrrus is directly determined using 
the low frequency data points and assuming a k~ 3 dependence. 

We clearly see in Figure 4 an excess over Pdrrus between fc=0.25 and 
0.6 arcmin -1 , which is more than a factor of 10 at k=0A arcmin -1 . Any reason- 
able power law spectrum for the cirrus component multiplied by the footprint 
leads, as can be easily seen in Fig. 4, to a very steep spectrum at spatial fre- 
quency k >0.2 arcmin -1 , which is very different from the observed spectrum. 
Moreover, the excess is more than 10 times larger than the measured instrumen- 
tal noise power spectrum. Therefore, as no other major source of fluctuations 
is expected at this wavelength, the large excess observed between k=0.25 and 



8 



Puget & Lagache 







, , | , , 


















-= 








E 


; 








r 














■ 


; 




















































***** 
































\ 

, , i , , * , 








0.1 








k (arcmin -1 ) 





Figure 4. Power spectrum of the source subtracted "Marano 1" field 
(*). The instrumental noise power spectrum (dotted line) has been 
subtracted. The dashed line represents the cirrus power spectrum mul- 
tiplied by the power spectrum of the PSF. 
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Figure 5. Power spectrum of the FIRBACK /ELAIS N2 field (o). 
The straight continuous and dashed lines, respectively, are the best 
fit cirrus power spectrum and the spectrum deduced from Miville- 
Deschenes et al (in preparation). The continuous curve shows the 
detector noise. 
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Figure 6. Power spectrum of the FIRBACK Nl extragalactic com- 
ponent. The continuous curve with o symbols represents the residual 
after the best-fit cirrus power spectrum has been subtracted from the 
power spectrum of the observational data in this field. The dotted 
lines represent the residuals after subtraction of the two power spectra 
of Miville-Deschenes et al., analogous to the two shown in Figure 5 
for the N2 field. The dot-dashed line is the detector noise spectrum. 
The dashed line is the residual extragalactic power spectrum obtained 
by subtracting 10 times the best-fit cirrus power spectrum from the 
directly observed power spectrum for this field. 



0.6 arcmin is interpreted as due to unresolved extragalactic sources. 

The Marano 1 field cannot be used to constrain the clustering of galaxies 
due to its rather small size. However, the extragalactic source power spectrum 
mean level can be determined. We obtain P SOU rces = 7400 Jy 2 /sr, which is in 
very good agreement with the spectrum predicted by Guiderdoni et al. (1997). 
This gives CFIRB rms fluctuations around 0.07 MJy/sr (for a range of spatial 
frequencies up to 5 arcmin -1 ). These fluctuations are at the ~9 percent level, 
which is very close to the predictions of Haiman & Knox (2000). 

The same analysis can be applied to the other and larger FIRBACK fields 
(Nl and N2). Fig. 5 shows the fluctuation power spectrum (P m ap) obtained for 
the FIRBACK /EL AIS N2 field in the plane of the sky. Again, we clearly see a 
large excess over P c %rrus- 

From Eq. 3, we deduce: 

Psources = {P-map Pnoise) /Wk ~ Pcirrus- 
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We see in the Nl and N2 fields (Fig. 6 and 7) an increase by a factor of 5 to 
10 of Psources from small to large scales. This increase cannot be due to cirrus: 
for example, the removal of 10 times the cirrus P/, leads to the same increase 
(dashed line in the two figures). So, it seems that we have with FIRBACK the 
first indication of IR sources clustering. Although the present analysis cannot 
be conclusive on this question, it clearly shows that the FIRBACK survey has 
the required sensitivity and is free enough from contamination to deal with this 
question. 

If we remove from the Nl and N2 maps all detected sources (S >200 mJy, 
5cr), we obtain CIB fluctuations which are compatible with a Poissonian distri- 
bution (Fig. 8). We need in fact in this case an independent cirrus template 
(like that for the HI gas) to remove the galactic contribution more correctly and 
see if correlations are present. 

In summary, we have seen that the fluctuations probe the CFIRB dominant 
source population. We detect the clustering of the detected FIRBACK sources. 
On the contrary, the dominant contribution to the CIB fluctuations seems to be 
compatible with a Poissonian distribution. 

It is important to note that the level of the extragalactic fluctuations de- 
rived from FIRBACK fields is in very good agreement with the Matsuhara et 
al. (2000, and this conference) study in the Lockman Hole region, where the 
interstellar dust contamination is negligible. The presence of correlated extra- 
galactic fluctuations should be studied using all the extragalactic ISOPHOT 
fields and independent cirrus tracers. 
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Figure 8. Power spectrum of the FIRBACK/ELAIS N2 unresolved 
extragalactic component (sources with flux S >200 mJy have been 
removed from the maps.) Line symbols are the same as in Fig. 6. 
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Discussion 

Martin Harwit: While the SCUBA sources might be predominantly at high 
redshifts, and cannot be ruled out as being at high z, the background is so low 
there that these galaxies are not major factors as far as production of heavy 
elements or massive star formation is concerned. Otherwise the abundance of 
chemical elements at high redshifts would be seen to be much higher. 

Jean-Loup Puget: Yes. In the model I presented the bulk of the energy 
comes from redshifts around z = 1, as do the heavy elements. 

Ray Norris: At higher redshift, metallicity will decrease, so presumably the 
dust fraction will decrease. To what extent can this account for your observed 
luminosity evolution? 
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Puget: The luminosity function I showed is a purely empirical one fitting 
both the background and the counts at 15, 170, and 850 fim. The dust contri- 
bution, which should decrease at large z, has to be combined with morphology 
of the galaxies. Local starbursts show much larger extinction than the average 
one. 



